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We describe a method to compute renormalized coupling constants in a Monte Carlo renornialization group 
calculation. It can be used, e.g., for lattice spin or gauge models. The basic idea is to simulate a joint system of 
block spins and canonical demons. Unlike the Microcanonical Renornialization Group of Creutz et al. it avoids 
systematical errors in small volumes. We present numerical results for the 0(3) nonlinear fr-model. 



In the Monte Carlo Renornialization Group 
(MCRG) [1] block spin transformations (EST) 
are analyzed by Monte Carlo (MC) simulations. 
In its original formulation this is done in terms 
of blocked observables suitable for the computa- 
tion of critical properties. This includes the de- 
termination of critical exponents from linearized 
renornialization group (RG) transformations [2] 
and the matching method for the calculation of 
the A/3 function [3]. Another formulation aims 
at the computation of effective actions generated 
by a EST. Here the difficulty is that effective ac- 
tions can be complicated with many couplings. In 
addition they can be hard to compute. 

Creutz et al. propose to compute renormalized 
couplings in a microcanonical demon simulation 
[4]. Their method rests on the equivalence of a 
microcanonical and a canonical ensemble in the 
thermodynamic limit. In small volumes it has 
systematic errors. This disadvantage can be over- 
come by making the demons canonical [7]. The 
resulting method requires also just a standard up- 
date program for its implementation. 

1. Description of the Method 

We consider a lattice spin system with spins 
(fix, action S, and Eoltzmann factor is exp(— 5). 
The action is assumed to have the form S = 
~'l2a l^a'^a with interactions Sa and couplings 
l3a. We then introduce an auxiliary demon sys- 
tem with action Sd = '^aPada and identical 



couplings Pa- The demons da are real numbers in 
the interval [0, dmax\- Consider the joint partition 
function 
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of spins and demons. It factorizes. The expecta- 
tion value of the demon variables d„ is 
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Once it is known we can compute /3„. 

Let us turn to a numerical RG transformation. 
We can simulate the effective theory without ex- 
plicit knowledge of the effective action. We simu- 
late a spin model with a given action, and apply 
a certain blocking rule to the configurations on 
the fine grid to produce block spin configurations 
on a coarser grid. The blocked configurations are 
distributed according to an effective Eoltzmann 
weight. The effective action can then be com- 
puted by a simulation of the above joint partition 
function. 

This joint simulation can be realized as follows. 
Assume that the effective action is of the above 
general form. First perform microcanonical up- 
dates of the joint system. These updates do not 
change the differences Sa — da- Since also the 
Eoltzmann weight is unchanged, knowledge of the 
l3a is not required for the update. Then replace 
the block spin configuration by a statistically in- 
dependent one. The second step is ergodic and 
fullfils detailed balance. Statistical independence 



of the block configurations can be assured if the 
configurations are sufficiently well separated in 
computer time. 

2. Numerical Results 

We implemented this method for the 0(3) in- 
variant vector model in two dimensions. We took 
an action S = — J2a=i Pa '^^ ^x,a with 12 inter- 
action terms S^^a = \ Y^y^Y^ „ {sx ■ Sy)" ■ The s^ 
are 3- vectors of unit length. Y^^a consist of all 
lattice points obtained by the obvious symmetry 
operations from a representative lattice vector v, 
and n takes values in {1, 2, 3}. Table 1 gives n and 
V for a = 1...12. As a first test we reproduced 

Table 1 

Labelling of the 12 interaction terms S^^a 
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given couplings. The model was simulated with 
some ad hoc chosen set of couplings on an 8^ and 
a 16^ lattice. The updates were performed with 
a tunable overrelaxation algorithm. The demon- 
spin update was done as follows. First, a new 
value for a single spin was proposed. Then we 
checked whether da plus the change of Sa was 
within [0,dmax]- If this was the case the system 
was updated. 

We employed 100 independent demon systems 
in order not to waste too many block spin config- 
urations. We generated a sequence of spin con- 
figurations, two successive configurations being 
separated by a number of sweeps large compared 
to the autocorrelation time divided by 100. The 
configurations for the 100 demon systems were 
then successively taken from this sequence. After 
a replacement of the spin configuration, we al- 
ways performed one lattice spin-demon updating 
sweep. After a full cycle through the demon sys- 
tems the measured demon values were averaged 
over the 100 systems and written to disk. 
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Reproduction of the original coupling constants 




couplings 


L = % 


i. = 16 


/3i 


1.30 


1.3010(23) 


1.2999(12) 


P2 


0.35 


0.3481(12) 


0.3485(8) 


P3 


0.01 


0.0105(9) 


0.0098(6) 


Pa 


0.02 


0.0190(7) 


0.0192(4) 


k 


0.004 


0.0045(6) 


0.0039(5) 


ii. 


-0.200 


-0.2045(30) 


-0.2002(18) 


P7 


-0.080 


-0.0808(21) 


-0.0808(12) 


Ps 


-0.020 


-0.0201(11) 


-0.0188(6) 


/39 


-0.01 


-0.0086(6) 


-0.0085(4) 


/3io 


-0.005 


-0.0045(8) 


-0.0038(6) 


/3ii 


0.02 


0.0248(26) 


0.0220(14) 


/3l2 


0.01 


0.0131(19) 


0.0114(12) 


Table 3 






Truncation experiment 






couplings 


L = % 


i. = 16 


Pi 


1.30 


1.1399(10) 


1.1408(7) 


P2 


0.35 


0.2993(6) 


0.3014(4) 


P3 


0.01 


- 0.0020(6) 


-0.0017(4) 


Pa 


0.02 


0.0138(3) 


0.0149(2) 


/35 


0.004 


0.0019(4) 


0.0022(2) 



Our results are presented in Table 2. In addi- 
tion to the reproduction of the 12 original cou- 
plings we studied the truncation to a subset of 5 
couplings (Table 3). For the L = 8 lattice we per- 
formed 20000 updates of the 100 demon systems. 
In the _L = 16 case we made 10000 updates. The 
results show the couplings are reproduced within 
error bars. 

For the original set of couplings we measured 
a correlation length ^ = 57.8(4). The correlation 
length of the truncated coupling set (L = 16) 
turned out to be smaller than the original one. 
We found ^trunc = 36.7(2). Note, however, that 
a naive truncation of the original set of couplings 
to the first five couplings leads to ^ = 195.5 ± 2.3 
on an _L = 400 squared lattice. 

Next we studied renormalization group trans- 
formations, starting from the standard action 
with nearest neighbour coupling only. We used a 
blocking rule that we call "dressed decimation" . 



Table 4 

Results for the first four RG steps, starting from the standard 

the correlation lengths for the effective theories. The columns g 



action with /3i = 
ive the results of 



= 1.9. The last line gives 
the subsequent RG steps 



RG step 





1 


2 


3 


4 


/3i 


1.9 


1.4619(15) 


1.1764(11) 


0.9785(16) 


0.8201(14) 


P2 








0.2691(9) 


0.3117(8) 


0.2892(8) 


0.2528(8) 


k 








-0.0123(6) 


0.0107(5) 


0.0226(6) 


0.0272(4) 


Pa 








0.0074(4) 


0.0132(3) 


0.0182(3) 


0.0195(4) 


k 








0.0043(5) 


0.0031(4) 


0.0032(5) 


0.0034(4) 


ii. 








-0.2213(19) 


-0.2611(13) 


-0.2257(14) 


-0.1689(12) 


P7 








-0.0800(12) 


-0.1027(10) 


-0.0985(13) 


-0.0766(9) 


Ps 








0.0003(9) 


-0.0057(6) 


-0.0104(7) 


-0.0107(8) 


Pg 








-0.0022(5) 


-0.0011(5) 


-0.0077(5) 


-0.0081(5) 


/3io 








-0.0005(6) 


-0.0011(5) 


-0.0014(7) 


-0.0015(7) 


/3ii 








0.0659(17) 


0.0898(13) 


0.0748(18) 


0.0470(15) 


/3l2 








0.0327(13) 


0.0440(11) 


0.0425(15) 


0.0301(12) 


e 


121.2(6) 


57.0(4) 


24.5(2) 


10.47(5) 


4.78(2) 



All lattice points that have coordinates x = (i,j) 
with i and j even are identified with block sites. 
The block spin s' at site x is then defined as 
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where the sum is over the nearest neighbours ofx. 
Tests with the massless free field theory revealed 
that the choice w = 0.8 is a good one in the sense 
that the effective actions and especially the fixed 
point action had good locality properties. 

As a first test we started with /3i = 1.9, and 
all the other couplings put to zero. The correla- 
tion length for this action is 121.2(6) [5]. We then 
blocked a 32^ lattice down to a 16^ lattice. The 
effective 12 couplings were then determined from 
the demon expectation values and used as input 
for the next iteration. Our results for the first 4 
steps are presented in Table 4. The last line shows 
correlation length estimates for the effective the- 
ories. The correlation length should theoretically 
be exactly halved each step. The significant devi- 
ations indicate that the action is not sufficiently 
well approximated by our ansatz. Note however 
that there is a very good decay of the couplings for 
given n with increasing distance of the spins. We 
conclude that local interactions with higher n and 
also interactions with more than two spins can- 
not be left out. (E.g., 4-point operators proved 



to be important in the study of Hasenfratz et al. 
[6] using the classical approximation.) 
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